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Until now, tlie generation of an excitonic reservoir in a cavity polariton system under quasi- 
resonant pumping lias always been neglected. We show that in microcavities having a small Rabi 
splitting (typically GaAs cavities with a single quantum well), this reservoir can be efficiently pop- 
ulated by polariton-phonon scattering. We consider the influence of the exciton reservoir on the 
energy shifts of the resonantly pumped polariton modes. We show that the presence of this reser- 
voir effectively reduces the spin anisotropy of the polariton-polariton interaction, in agreement with 
recent experimental measurements, where the multistability of cavity polaritons has been analyzed 
[Nature Materials 9, 655-660 (2010)]. 
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I. INTRODUCTION 

Exciton-polaritons (polaritons) are the quasiparticlcs 
formed of photons and excitons strongly coupled in 
microcavitiesi. Their popularity comes from their unique 
properties, mixing those of light (small effective mass) 
and excitations of matter (self-interactions and thermal- 
ization with phonons). Polaritons are photonic parti- 
cles interacting with each other via a well-defined mi- 
croscopic mechanism: exciton-exciton interaction. As a 
result, strongly coupled microcavities show a very strong 
non-linear optical response. Parametric oscillations and 
amplification^rii have been demonstrated a decade ago, 
showing a record efficiency, whereas low-threshold bista- 
bility has been demonstrated in 2004^. Polaritons also 
possess two spin projections which can lead to non- 
linear spin dependent effects, allowing to use microcav- 
ities as polarization-selective optoelectronic devices, or 
" spin-optronic" devices^. The first method of controlling 
the polariton spin degree of freedom is the splitting exist- 
ing between the TE and TM polarized polariton modes. 
The effect of this splitting can be mathematically de- 
scribed as an effective magnetic field acting on the po- 
lariton pseudospin. It is at the origin of time-dependent 
spin oscillations^ and of the optical spin Hall effect^i^. 
The spin degree of freedom can be used to generate half- 
integer topological defects such as half-solitonsi^. 

Another key ingredient of the polariton spin dynamics 
is the spin anisotropy of the polariton-polariton interac- 
tion which was noticed for the first time ir>ii. Exciton- 
exciton interaction constant in the triplet configuration 
ai is dominated by Coulomb exchange interaction. It 
can be calculated in the Born approximation similarly 
to the spin-less caso^. On the other hand, an exchange 
process between excitons having opposite spins ±1 (sin- 
glet configuration) leads to the formation of dark exci- 
tons having a spin ±2. The polariton-polariton interac- 
tion in the singlet configuration is therefore a second- 
order process which invokes intermediate states (typi- 



cally located half a Rabi splitting above the polariton 
state). The resulting interaction is weakly attractive 
(constant 02), much weaker than the interaction in the 
triplet configurationi^ii^. This spin anisotropy led to the 
prediction and observation of a wide variety of remark- 
able phenomena, such as the inversion of the linear polar- 
ization during polariton-polariton scattering processes^^, 
the spin-Meissner effeciii^, the stability of half-topological 
defectsi^, and the polariton multistabilityiiii^, which has 
been recently evidenced experimentallj*i^"— . The sur- 
prising conclusion of the latter experiments^ was that the 
theoretical description based on the driven-dissipative 
Gross-Pitaevskii equation required to consider a repulsive 
polariton-polariton interaction in the singlet configura- 
tion in order to fit the experiments. This conclusion was 
also based on recent observations suggesting that the bi- 
exciton state could play the role of the intermediate state 
for the second order polariton-polariton interaction^^. As 
a result, a2 could become comparable with ai and change 
sign when polariton energy becomes larger than the bi- 
exciton energy. 

It is worth to mention here that a description of 
the resonantly pumped polariton system by the Gross- 
Pitaevskii equation, widely used nowadays, is based on 
the assumption that the majority of the polaritons occu- 
pies the coherent states. However, the systematic time- 
resolved transmission experiments in ns time domai n^^i^" 
have clearly shown the important role of long-living 
states (~ 500 ps lifetime) populated under resonant ex- 
citation. The experimentally observed blue shift of the 
transmitted light was shown to be different from the tem- 
poral profile of coherent polariton density. The blue shift 
remained of the order of meV well after the coherent 
polaritons have decayed. This type of effect cannot be 
described within the Gross-Pitaevskii approach, where 
the only source of energy shifts are the coherent polari- 
tons. In order to explain these experimental findings, the 
authors have considered phenomenological coupling of 
the resonantly pumped polariton mode with a long-lived 
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"reservoir" states. Although the exact nature of both 
" reservoir" and coupHng was not specified, the success in 
the description of the wide experimental data set (also for 
the case of cw resonant pumping2i) has proven the im- 
portant role of the "reservoir" in these experiments. This 
concept also allowed to resolve several contradictions in 
the interpretation of data obtained in cw and ps-range 
pulsed excitation experiments (see detailed discussion in 
Refi^). A coupling with the reservoir has also had to 
be introduced to describe OPO polarization dynamics in 
earlier experiments^ 
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FIG. 1: (Color online) The scheme of the process we con- 
sidered. The circularly polarized pumping is in the ground 
state. The interaction with phonons provides an exchange of 
polaritons between the condensate and the reservoir. 

In the present work, we propose to describe theoret- 
ically the thermal generation of an excitonic reservoir 
from a resonantly pumped low energy polariton state. To 
our knowledge, this effect has never been studied, being 
always considered as negligible. As we are going to show 
below, this is essentially incorrect in structures showing a 
small Rabi splitting at zero or positive detunings, namely 
in the widely-studied GaAs-based cavities containing a 
single QW. This is the type of structure which has been 
used so far in multistability experiments^ii^. In the 
first part of the paper we solve the semi-classical Boltz- 
mann equations with the ground state being continuously 
pumped by a resonant, circularly polarized (cr^) laser as 
sketched on the figure 1. This way of using Boltzmann 



equations is extremely unusual for the field of exciton- 
polaritons, since such equations are typically used in or- 
der to describe the non-resonant pumping case and the 
relaxation of reservoir excitons toward the lower energy 
states^. We show that in the above-mentioned type of 
structures, the resonant pumping results in the forma- 
tion of a thermal exciton reservoir of moderate density 
which we assume to be unpolarized whatever the polar- 
ization of the pumping laser. This assumption might not 
be fully correct, since the interactions with the strongly 
polarized pumped should split the circularly polarized 
reservoir state, which might slow down the standard ex- 
citon spin relaxation mechanisms^^. In the absence of 
any data allowing to verify the presence of such effect, 
we are going to assume throughout this paper that the 
exciton spin relaxation is much faster than the reser- 
voir exciton lifetime, leading to an unpolarised reservoir. 
This reservoir can provoke a blue shift of both polari- 
ton spin components by a fraction of meV, comparable 
with the energy shifts induced by the resonantly pumped 
state occupation. This can result in an apparent posi- 
tive a2, induced by the reservoir. In the second part of 
the paper we solve sclf-consistently the coupled spinor 
Gross-Pitaevskii-Boltzmann equations, using the proper 
attractive 012 constant for polariton-polariton interaction. 
We show that the correct accounting for the reservoir 
generation by acoustic phonons is sufficient to describe 
the most important features of the experiment of Refi^. 

II. ANALYTICAL APPROXIMATION 

We begin with a simple analytical model allowing 
to estimate the importance of the effect we would like 
to discuss, that is, the population of an exciton reser- 
voir by scattering with phonons up from the coherently 
pumped condensate. The spatially homogeneous conden- 
sate and the reservoir are described by the occupation 
numbers. In what follows we will use lowercase letters to 
describe the surface densities (e.g. the condensate den- 
sity no = Nq/S) and uppercase ones for the occupancy 
of the states. 

In the most general case the semi-classical Boltzmann 
equation for the occupation number of the state with 
wave-vector k reads^^: 



dt 



k' 



(1) 



Here the terms with W are the scattering rates between 
different states, P is the pumping and T are the decay 
rates of the corresponding states. The terms (I + N£) 



account for the bosonic nature of the particles involved, 
allowing to describe bosonic stimulation towards strongly 
populated states. 
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This system of equations will be solved as a whole be- 
low, but first we will discuss an approximated solution 
allowing analytical treatment. Indeed, if one assumes 
the conditions when single-phonon assisted scattering 
is stronger than multi-scattering channels, the reservoir 
states become only coupled with the ground state, and 
not between each other. In QW-based structures, the 
exciton-phonon coupling arises between 2D excitons and 
3D phonons. The wave vector in the plane is there- 
fore conserved, which is not the case for the wave vector 
of phonons in the z direction (g^). This specificity al- 
lows to organize energy-conserving scattering processes 
even in the very sharp polariton dispersion. This process 



is however limited by a wave vector cutoff provided by 
qz,max ~ L being the well width. In practice this 
cutoff limits the energy exchanged during a polariton- 
acoustic phonon scattering event to 2-3 meV. So qualita- 
tively, one expects single-scattering processes to be dom- 
inant if the energy difference between polariton ground 
state and the excitonic reservoir is smaller than 3 meV, 
and multi-scattering processes if this energy difference is 
larger than 3 meV. 



In the framework of the above-mentioned assumptions, 
the equation of motion for the reservoir state k writes: 



J 



dt 



-r^iV,- - Nf:Wj:^„{l + N„) + (1 + Ni:)W,^j:No 



(2) 



Using the Fermi golden rule and assuming Lorentzian 
broadening of polariton lines, one can calculate the scat- 



tering rate for polariton-acoustic phonon interaction: 
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(3) 



Here 7^ is the broadening of the polariton line corre- 
sponding to the decay rate in the Boltzmann equa- 



tions, LOq - the phonon frequency, N 



ph 

ber of phonons, M{q) - the matrix element of polariton- 
phonon interaction, taking into account the excitonic 
fraction of each of the states involved. In order to fur- 
ther simplify this equation, we assume that the only 
dependence of WS'^°p between the ground state and 

the reservoir state is via the number of phonons N!^ ■ 



ph ■ 



WN^^ and Wj:_ 



1) where N!^ 



0-^k " "'ph """" " k^O " y^'ph ~^ ) "^i^i-- ^i ph 

l/{exp{-AEk/kBT) - 1). We take the value of W as 
a constant, which will be shown in section IV to pos- 
sess a value of order 10^s~^ using an exciton decay rate 
Ffc w h/AOQ ps. Thus, the above-mentioned wavevec- 
tor cutoff is neglected together with the wavevector de- 
pendence of the interaction with acoustic phonons, and 
the only dependence that is kept is the exponential de- 
crease of phonon occupation numbers for large energies 
exchanged. 



For this system of equations, one can analytically ob- 
tain the stationary values of the occupation numbers: 



N^h 



Nn 



(4) 



from which the total reservoir density can be straight- 



W 



critical 



forwardly obtained. Here N§ = 
value for No, discussed below. 

From the previous equation, one can see the existence 
of two distinct regimes. If Nq is smaller than N§, the 
reservoir density increases linearly with A^o- It then sat- 
urates when Nq ^ Nq. In the latter regime, the reser- 
voir density does not depend on pumping anymore and 
is proportional to the phonon distribution function. In 
this regime, the reservoir distribution function is there- 
fore exactly by a Bose distribution function with chem- 
ical potential equal 0, if the energy is counted from the 
resonantly pumped polariton mode. 

The results of the calculations are presented on the 
figure 2. The dashed lines show the reservoir density 
calculated with the procedure described above in the sat- 
urated regime versus the exciton-photon detuning for 3 
different temperatures. The Rabi splitting of the struc- 
ture is 4 meV, which corresponds to the case of a GaAs 
cavity containing a single QW— Going towards pos- 
itive detuning reduces the energy splitting between the 
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ground state and the reservoir, which increases the num- 
ber of phonons. One can also increase the number of 
phonons directly, by simply increasing the temperature. 
In both cases, the reservoir density increases. 

Next, in order to check the validity range of our analyt- 
ical approximation, we have solved numerically the whole 
system of semi-classical Boltzmann equations where all 
possible scattering paths are taken into account Eq|T] 
The pumping term is different from zero only for ground 
state Pq 7^ 0- The simulation is carried out until the 
system reaches a steady-state distribution Nj^. 

The results of the numerical simulations are presented 
on the figure 2 with solid curves. As expected, the analyt- 
ical approach shows a good accuracy in the positive de- 
tuning range when the energy difference between the po- 
lariton ground state and the excitonic reservoir is small. 
In that range, and at 20 K, the densities found arc likely 
to destroy the strong coupling regime, and anyway to pro- 
voke strong blue shift due to exciton-exciton interactions. 
By going to negative detuning, the generated reservoir 
density drops very rapidly, becoming completely negligi- 
ble. This drop is explained by the fast decay of the direct 
scattering processes and the inefficiency of the multiple 
scattering channels. One should also note that the ac- 
counting for polariton-polariton scattering mechanisms 
brings no change to the present picture as soon as only 
the polariton ground state is pumped. 
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FIG. 2: (Color online) Densities of excitons versus detuning 
for different values of temperature, solid lines - numerically 
calculated, dashed lines - analytical approximation. Open 
circles - density of condensate no. Rabi splitting 4 meV. 

The simulations using the complete system of Boltz- 
mann equations have also been performed for a GaAs 
microcavity showing a significantly larger Rabi splitting 
of 12 meV. The results are presented on the Fig. [31 In 
that case, the reservoir density generated by the resonant 
pumping is completely negligible, keeping below to 10^ 
cm~^ for the values of detuning considered, and the an- 
alytical approximation of singlc-phonon scattering is not 
applicable. 
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FIG. 3: (Color online) Densities of excitons versus detuning 
for different values of temperature. Rabi splitting 12 meV 



III. INFLUENCE OF THE RESERVOIR ON 
ENERGY SHIFTING 

Our next step is now to consider the impact of the ex- 
citonic reservoir, generated under quasi-resonant pump- 
ing, on the energy of the two spin components of the 
condensate. The interaction between excitons has been 
previously studied theoretically and experimentally, both 
in the scalar— and spinor cases^^^, taking the account 
the fact that the excitons are composite bosonsi^. We 
assume that this interaction is the same, whether the ex- 
citons concerned are both in the condensate, or one of 
them is in the reservoir. 

In the triplet configuration the 2D interaction constant 
comes mostly from the exchange Coulomb interaction. It 
can be calculated in the Born approximation yielding the 
approximate formula ai = 3Ei,a^. In the singlet config- 
uration, the exciton-exciton interaction is a second-order 
process passing either through the dark exciton states 
possessing spin ±2 or trough bi-exciton states. These 
latter states are usually several meV above the polari- 
ton states, which has two consequences: the matrix ele- 
ment of interaction becomes reduced and obtains a neg- 
ative sign, as any second-order correction for the ground 
state. We should note, however, that it is possible to have 
the biexciton state close to the polariton ground state, if 
the Rabi splitting is small and the detuning is strongly 
positive^. In what follows, we assume that the singlet 
interaction constant is a2 = — O.lai. 

The corrections of the energies of the two polarization 
components of the condensate take into account both in- 
teraction constants: 

AE'^ = aiTi^ + Q!2?T._ (5) 
AE~ = airt+ — a2n^ (6) 

Where £'+'^ arc the energy shifts of the (7+ and (t~ 
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polarized component, are the exciton densities (in- 

cluding both the condensate and the reservoir) and ai^2, 
the interaction constants. If we decompose the total exci- 
ton density into reservoir density (the reservoir is not 
polarized) and condensate density no (the condensate is 
polarized ct^, as the pump), we can write the energy 
shifts as follows: 

= a2y +ai (y +no) (7) 

A^^ =a2 (y+Jloj +aiy (8) 

Using this formula, we plot on the Figure 4 the en- 
ergy shifts of the two components of the condensate as a 
function of the detuning, using the values of the reservoir 
density obtained by solving the full Boltzmann equations. 
The results presented on the figure correspond to a Rabi 
splitting of 4 meV and a temperature 5K. 
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FIG. 4: (Color online) Energy shift of cr+ (red) and cr~ (black) 
components of condensate . Rabi splitting 4 meV, tempera- 
ture 5K. 

We see that in the conditions of this simulation, 
namely, at relatively small Rabi splittings and quasi- 
resonant cw pumping (as opposed to pulsed excitation), 
the contribution of the reservoir can totally modify the 
observed energy shifts. Indeed, at negative detuning the 
reservoir is not strongly populated, and a redshift of the 
a~ component is observed, as one would naturally expect 



from the negative sign of the a2 constant. However, at 
positive detuning the reservoir starts to play a role. Since 
the reservoir is depolarized, it contains a large number of 
a~ excitons which induce a blue shift for the cr" compo- 
nent of the condensate via the stronger repulsive interac- 
tions, which dominate the weak redshift induced by the 
excitons. 

Therefore, a measurement of the energy shift of the 
component at positive detunings, without including the 
filling of long-lived reservoir states, can lead to inexact 
conclusions concerning the wrong sign and the magnitude 
of the polariton-polariton interaction constants. 

IV. APPLICATION OF THE MODEL 

One of the important effects based on the energy shifts 
of the polariton mode is the bistability (or multistability 
in the spinor case) of the polariton mode under quasi- 
resonant pumping. Indeed, if the system is pumped with 
a cw laser whose frequency is above that of a bare po- 
lariton mode, blue shift of the macrooccupied mode will 
bring it closer to the laser, increasing the efficiency of 
the pumping. Therefore, at some pumping intensity the 
mode becomes unstable and the population jumps up 
abruptly. This effect has been predicted to occur also 
in the spinor case^, but if one takes into account the 
negative sign of 012, it is logical to expect that the jump 
of one polarization component under elliptical pumping 
will prevent the jump of the other component by moving 
it off-resonance. 

In the first experimental work on polariton 
multistability^, the authors have, in particular, 
observed an abrupt simultaneous increase of the ampli- 
tudes of both components, when the polarization of the 
excitation light is almost linear, but slightly elliptic (Fig 
l(d,e) in Ref.— ). This contradicts the expected behavior 
mentioned above, and was explained by the authors 
by introducing a positive sign of 02, that is, repulsive 
interactions for excitons with opposite spins. 

In this section we show how the formation of an exci- 
tonic reservoir in this particular type of experiment in- 
fiuences the observed behavior of the polariton modes in 
the multistable regime. In order to describe correctly 
the behavior of the macrooccupied modes under quasi- 
resonant pumping, we write the coupled Gross-Pitaevskii 
equations with pumping and decay for both polariza- 
tion components taking into account the interactions 
with the reservoir. 
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(9) 



Here, 7 defines the broadening of the polariton line, S density and n^, the density of excitons in the reservoir, 
the surface area of the polariton spot, P±, the excitation The detuning between the laser and the pumped polari- 
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ton state A = 0.5meV, as in the experiment. The density 
of the reservoir can be obtained analytically, following the 
approach introduced in section I: 



1 

2^ 



iV^\(|*+|2 + |vI,_|2) 



1 



■kdk (10) 



Taking an average value for and considering 

parabolic exciton dispersion, we can evaluate this integral 
analytically and obtain the solution for function 
of (|*+|2 + |vl/_|2): 



A = 



+ If^P + 1 



B 



(11) 
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By introducing the above expression in the eq. [9l we 
get a closed equation for the amplitude of the pumped 
polariton mode. An essential parameter of the model is 
the ratio (which has to be assumed constant in or- 
der to obtain analytical results). It can be estimated by 
calculating with the full system of semi-classical Boltz- 
mann equations the value of the reservoir density Ur as 
a function of the condensate density no — i£±L^i£.^^ 
We have then fitted the results obtained numerically by 
the analytical formula [TT] The results of this fitting are 
presented on Fig. 5. The value obtained for the fitting 
parameter was around 300. 



By changing the pump density slowly in time from to 
some value and back again we can obtain the evolution 
of the populations of both circularly polarized polariton 
components and of the excitonic reservoir. The results 
are presented on the Fig. 6. Indeed, a simultaneous jump 
of both circular components, and of the reservoir is visi- 
ble when increasing the pumping power. The strong rise 
of the cr"^ polariton density associated with the passage 
of the bistable threshold leads to a strong rise of the 
population of the unpolarized excitonic reservoir which 
provokes a blue shift of the polariton energy of the a~ 
component. This blue shift overcomes by far the red shift 
induced by the attractive interaction between the oppo- 
site spin components of the condensate. It leads to the 
passage of the bistable threshold for the a~ component 
as well. 
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FIG. 6: (Color online) Dependence of the amplitudes of 
(red) and (black) components on pump power. Polariza- 



tion of pumping pp^ 
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FIG. 5: (Color online) Black circles - dependence of Ur on 
no calculated numerically. Red curve - fit of this data by the 
eqs. 111! Fitting gives = 274. This value corresponds to 

"-8 = ^ = 3- 10^ 



Another experimental effect is well reproduced when 
the system is on the upper stability branch and the 
pumping power is decreased. The cr" -component has a 
smaller population and therefore jumps down first. At 
the same time, the population of the other component 
strongly increases. This effect is easy to understand 
as well. When a~ jumps down, the energy of cr+ is still 
higher that the laser frequency. The jump of leads 
to a reduction of the total reservoir density n^, bringing 
the frequency of the a'^ component closer to the laser 
frequency. This makes the pumping more efficient and 
increases the (7+ population. When this amplitude in 
its turn abruptly decreases, the reservoir density jumps 
down once more, and this causes a second jump down of 
the (T~ -component. So, while the pumping power is de- 
creased, the amplitude of (T~-component undergoes two 
steps, as it was observed in experiment. 

However, in the framework of this model, we have 
not been able to fully describe the experimental behav- 
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ior reported ins^ for large circular polarization degree 
of the pump. We believe that the hypothesis of the 
fully unpolarized excitonic reservoir is no more valid in 
the case of quasi-circular pumping. Indeed, the jump 
up of one circular component only induces an energy 
splitting between and ct+ excitons in the reservoir, 
which may substantially slow down the spin relaxation. 
Thus, the reservoir may become spin polarized which 
would strongly affect the behavior of the pumped modes. 
The proper analysis of these phenomena requires a more 
complete description of the dynamics of the coupled 
condensate-reservoir system and will be addressed in fu- 
ture works. 

The authors acknowledge the support of EU ITNs 
"Spin-Optronics" 237252 and "INDEX" 289968, IRSES 
"Polaphcn" 246912 and ANR "Quandyde". We would 
like to thank Dr. Mikhail Glazov for useful discussions. 



of a polariton mode. Exciton-polaritons from the con- 
densate are scattered up with phonons. It was shown 
that in structures with small Rabi splitting and pos- 
itive exciton-photon detuning, this process provides a 
strongly populated reservoir which influences the behav- 
ior of the resonantly pumped polariton modes. Because 
of the reservoir, it is possible to observe a blueshift of 
the polariton condensate component polarized opposite 
to excitation, while the exciton-exciton interaction con- 
stant for opposite spins remains negative (attractive in- 
teraction). An analytical model allowed us to describe 
the multistable behavior of the polariton system recently 
evidenced experimentall y^ \^22 ^ 



V. CONCLUSION 



In this work we have studied the generation of an ex- 
citonic reservoir in the case of quasi-resonant pumping 
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